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^ (57) Abstract: An apparatus and method of suppressing seismic multiples (figure 1), used in the collection of marine seismic data 
2 ( 106 > witn tne sources (104), and receivers (105), positioned in the aquatic marine seismic survey environment (102), and with the 

receiver positioned in the water and near the below surface (103), and sources with a bubble diffuser (101), positioned so that the 
Q bubbles emitted by the diffuser are positioned between the receivers and the air-water surface, and wherein the emitted bubbles 

provide high acoustic reflection and suppress specular reflection of seismic waves by using a bubble control (106), to activate the 
^ diffuser during the collection of seismic data during the survey. 
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METHOD FOR MULTIPLE SUPPRESSION 
BASED ON PHASE ARRAYS 



FIELD OF THE INVENTION 

This invention relates generally to the field of geophysical prospecting. More 
particularly, the invention is a method of suppressing multiples in the collection of 
marine seismic data. 

BACKGROUND OF THE INVENTION 

In marine seismic data, one of the most significant sources of noise is 
reverberation caused by up-coming seismic waves that have previously reflected off 
the air-water interface. This type of noise is called 'multiples 1 and it re-illuminates 
water-bottom and sub-surface reflectors with time latency. Consequently, multiples 
generate false images of the sub-surface and hamper accurate interpretation of marine 
seismic data. The purpose of the present invention is to suppress or reduce multiple 
energy at the time of data collection 

There exist several data processing methods to eliminate or reduce the 
multiples in seismic data. Most of these methods are based on filtering multiples out 
of the data by taking advantage of their unique properties. These properties include 
predictive repetition with a constant amplitude decay rate and move-out difference 
from the primary signal. Other more advanced data processing methods involve 
modeling the multiples and solving the full wave equation. 

Another class of methods to reduce multiples in intermediate depths of water 
uses a configuration of marine seismic data acquisition called 'ocean bottom cable 1 
(OBC). This type of acquisition positions different sensors, typically hydrophones 
and geophones, in a cable at the water bottom to measure both pressure and particle 
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velocity waves, respectively. Combined with a proper processing algorithm, these 
methods can subtract a portion of multiple energy from the collected seismic data. 

Multiple contamination is a fundamental impediment to use of detailed analysis of 
seismic amplitude data in many prospect areas. Complex multiples severely corrupt 
target horizons rendering structural and basic stratigraphic mapping difficult and 
seismic amplitudes of littie value. Aggressive application of current mitigation 
methods has failed to substantively improve data quality, especially when the 
multiples are more than an order of magnitude stronger than the primary signals of 
interest. 

W. H. Clark's U.S. Patent No. 4,625,302, "Acoustic Lens for Marine Seismic 
Data Multiple Reflection Noise Reduction", issued Nov. 25, 1986, discloses a method 
for multiple reduction in marine seismic data acquisition using air bubbles. Clark's 
'302 patent relies on a refraction mechanism achieved by lowering local sound speed 
relative to nominal water sound speed. This is accomplished by a bubble emitter 
configured to produce an acoustic lens of air bubbles. The air bubble lens is designed 
to refract, rather than reflect, the multiples laterally and thus away from the receivers. 

P. D'Antonio et al. in U.S. Patent No. 5,401,921 "Two-Dimensional Primitive Root 
Difiusor 11 , issued Mar. 28, 1995 to RPG Difiusor Systems, Inc., discloses a two- 
dimensional acoustic diffuser based on a primitive root sequence. This diffuser is 
constructed with solid materials, not air-bubbles. It has a relatively small scale (2ft by 
2ft for one period of two-dimensional primitive root sequence) because it is designed 
for rfiffiising sound wavelengths in air. In addition to the phase grating, it also relies 
on sound absorption by the back plane. The main application is for suppressing sound 
in confined air spaces such as music halls and recording studios. 

W. H. Berhens' U.S. Patent No. 5,959,938, "Tuned Bubble Attenuator for 
Towed Seismic Array", issued Sep. 28, 1999, discloses a method to attenuate 
horizontally transmitted acoustic energy in a marine environment by emitted bubbles. 
Berhens' '938 patent uses the attenuation mode, rather than the reflection mode, by 
attempting to create relatively large Q> 9 mm) bubbles at the resonant frequency of the 
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seismic wave. The emitted bubbles are localized near the seismic source and their 
purpose is to restrict side-directional propagation to enhance signal-to-noise ratio, 
rather than to attenuate multiples. 

{SUMMARY OF THE INVEfrlTION 

The present invention is an apparatus for suppressing multiples in the 
collection of marine seismic data. The apparatus comprises at least one source 
positioned in the body of water, at least one receiver positioned in the body of water 
below the air-water interface and near the sources; a bubble diffuser positioned in the 
body of water so that the bubbles emitted by the bubble diffuser are positioned 
between the receivers and the air-water interface, wherein the emitted bubbles provide 
high acoustic reflection and substantially suppress specular reflection of seismic 
waves; and a control for activating the bubble diffuser during the collection of seismic 
data by the sources and receivers. 

The present invention is also a method for suppressing multiples in the 
collection of marine seismic data. At least one source and at least one receiver are 
positioned in a body of water, with the receivers positioned below the air-water 
interface. A bubble difiuser is positioned in the body of water so that the bubbles 
emitted by the bubble diffuser are positioned between the receivers and the air-water 
interface. The emitted bubbles provide high acoustic reflection and substantially 
suppress specular reflection of seismic waves. The bubble diffuser is activated during 
the collection of seismic data by the sources and receivers. 

BRIEF DESCRI PTION OF THE DRAWINGS 

The present invention and its advantages may be more easily understood 
by reference to the following detailed description and the attached drawings in which: 

FIGS, la and lb show top and side views, respectively, of an embodiment 
of the apparatus of the present invention for suppressing multiples in the collection of 
marine seismic data; 



WO 02/50755 



-4- 



PCT7US01/47837 



FIGS. 2a through 2d show top and stern views of two embodiments of the 
apparatus of the present invention; 

FIG. 3a shows an embodiment of the wedge shaped bubble diffuser with 
eight bundles of porous air hoses stacked vertically; 

FIG. 3b shows an embodiment of the wedge shaped bubble diffuser with 
ten air hoses wrapped around a supporting steel cable; 

FIG. 3c shows an embodiment of the wedge shaped bubble diffuser with 
300 trailing air lines; 

FIG. 4 is a flow chart illustrating the steps of the method for one 
embodiment of the present invention for suppressing multiples in the collection of 
marine seismic data; 

FIG. 5 illustrates an example of a vertically oriented phase grating 
constructed in accordance with the present invention; 

FIGS. 6a and 6b show profiles of a phase array computed for a primitive 
root sequence without and with, respectively, the focusing technique of the present 
invention; 

FIG. 7 shows the normalized amplitude of the transfer function of specular 
scattering for the multiple suppresser, 

FIG. 8 shows the comparison between the spectra of specular reflection for 
the input pulse without and with, respectively, the multiple suppresser; 

FIG. 9 shows the acoustic responses in the time domain of the multiple 
without and with, respectively, the multiple suppresser; 

FIG. 10 shows a wedge shaped embodiment of the multiple suppresser 
used in the small-scale field experiment example; 

FIG. 11 shows a primitive root embodiment of the multiple suppresser 
used in the small-scale field experiment example; 

FIGS. 12a and 12b show seismic images taken without and with, 
respectively, the wedge shaped multiple suppresser of FIG. 10; and 
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FIGS. 13a and 13b show seismic images taken without and with, 
respectively, the primitive root multiple suppressor of FIG. 1 1 . 

While the invention will be described in connection with its preferred 
embodiments, it will be understood that the invention is not limited thereto. On the 
contrary, it is intended to cover all alternatives, modifications and equivalents that 
may be included within the spirit and scope of the invention, as defined by the 
appended claims. 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention is an apparatus and method for suppressing multiples in 
the collection of marine seismic data. Unlike the previously mentioned data 
processing methods, the present invention directly attacks the generation of multiples 
at the air-water interface. The present invention does this by creating phase changes 
in the down-going waveform at the air-water interface. The physical concept of a 
phase array is known in the fields of antenna theory, radar, and sonar. By illu m i n a tin g 
an appropriately designed phase array, patterns of electromagnetic or acoustical fields 
can be manipulated. For example, phase arrays have been used to enhance forward- 
to-backward gain or angular sensitivity. Depending on the choice of phase grating, it 
is possible to reduce the specular reflection of the wave field and divert wave energy 
to directions of no interest. This type of phase grating is beneficial to suppress the 
multiple, which is generated by the specular reflection of up-going seismic waves at 
the air-water interface. 

The multiple-suppressing phase array of the present invention is physically 
implemented by creating surfaces of high acoustic reflection with depth variation near 
the air-water interface. Considering the range of wavelengths of interest for marine 
seismic data acquisition, it is not practical to use a solid structure with the required 
acoustical properties (e.g. near 100% reflection) to build a phase array. However, it is 
possible to reflect seismic waves to a desirable direction by injecting air bubbles into 
water. The present invention uses air bubbles (or a mixture of air and water) to 
achieve the necessary acoustic properties for the reflection of the seismic wave. This 
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is discussed in Domenico, S. N., "Acoustic wave propagation in air-bubble curtain in 
water - part I: History and theory", Geophysics, v.47, pp.345-353, 1982. A small 
amount of air in the water decreases sound wave speed so much that the seismic wave 
is reflected from the boundary of air-mixed water. The present invention can be 
applied to various configurations of marine seismic data acquisition including, but not 
limited to, towed single/multi streamer and 2 or 4 component ocean bottom 
cables/sensors (2/4C OBC/OBS). 

FIGS, la and lb show top and side views, respectively, of an embodiment of 
the apparatus of the present invention for suppressing multiples in the collection of 
marine seismic data using a specially configured bubble difiiiser 101. Specific 
embodiments for special configurations of the bubble difiiiser will be described in 
detail in conjunction with the discussion of FIGS. 2c and 2d, below. The bubble 
difiiiser 101 is positioned in a body of water 102 with an air-water interface 103. At 
least one seismic source 104 is positioned in the body of water 102 at or below the air- 
water interface 103, as is typical in the industry, and in the general vicinity of the 
bubble difiiiser 101. The seismic sources 104 are preferably airguns or marine 
vibrators, positioned singly or in arrays, but may be any other appropriate marine 
seismic source known in the art The seismic sources 104 will be used in conjunction 
with the bubble difiiiser 101 in the collection of marine seismic data. 

At least one seismic receiver 105 is positioned in the body of water 102 below 
the air-water interface 103 and in the general vicinity of the sources 104, as is typical 
in the industry. The seismic receivers 105 are preferably pressure sensors such as 
hydrophones, but may be any other appropriate sensors adapted for the marine 
environment. The seismic receivers 105 will be used in conjunction with the sources 
104 and the bubble difiiiser 101 in the collection of marine seismic data. 

In a preferred embodiment, both the sources 104 and receivers 105 are towed 
through the body of water 102 by a towing vessel 106. This towing vessel 106 is 
preferably a seismic vessel equipped for this task. The sources 104 and receivers 105 
may be towed by the same or different towing vessels 106. The sources 104 may be 
towed in single or multiple arrays and the receivers 105 may be towed in single or 



WO 02/50755 PCT/US01/47837 

-7- 

multiple seismic streamers. In an alternative embodiment, the receivers 105 may be 
positioned on the bottom of the body of water 102 in an ocean bottom cable. Here the 
receivers 105 are preferably pressure and particle velocity sensors such as 
hydrophones and geophones, respectively, but may be any other appropriate sensors, 
such as accelerometers, adopted for the marine environment. 

The bubble diffuser 101 is positioned in the body of water 102 at or below the 
air-water interface 103 and in the general vicinity of the sources 104 and the receivers 
105. The bubble diffuser 101 may be attached to a support 107 in some embodiments. 
The bubble diffuser 101 emits bubbles 108 into the body of water 102. The bubble 
diffuser 101 is positioned so that the emitted bubbles 108 are positioned between the 
receivers 105 and the air-water interface 103. 

FIGS. 2a through 2d show top and stern views of two embodiments of the 
apparatus of the present invention. A towing line 109 is attached to the bubble diffuser 
101 or the support 107 for the bubble diffuser 101 and then to a vessel for towing. 
The towing line 109 may be attached to and towed by the same towing vessels 106 
towing the sources 104 or receivers 105. Alternatively, the towing line 109 may be 
attached to and towed by auxiliary work vessels 110, such as one or more tugs, as 
shown in FIGS. 2a and 2b. The auxiliary work vessels 1 10 assist the towing vessels 
106 to deploy, tow, and recover the bubble diffuser 101. Alternative embodiments are 
also possible. For example, the bubble diffuser 101 may be towed in front of the 
towing vessels 106, rather than between the seismic sources 104 and receivers 105. 

A compressed air hose (not shown) for supplying compressed air to the bubble 
diffuser 101 is attached to the bubble diffuser 101 and to a source of compressed air 
(not shown). The source of compressed air may be located on the towing vessel 106 
or auxiliary work vessel 1 10 towing the bubble diffuser 101, but may also be located 
on a separate vessel (not shown). The compressed air hose is preferably designed to 
supply sufficient pressure so that the emitted bubbles 108 are of a sufficiently high 
volume concentration in the water to act as a high acoustic reflector of seismic waves. 
The bubble diffuser 101 is also preferably designed so that the emitted bubbles 108 
are of sufficiently small size that the bubbles 108 rise slowly in the body of water 102. 
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There are many ways to mix air into water to generate bubbles. For one 
example, pressurized air can be fed to a porous ceramic stone to generate fine air 
bubbles in the water. Another example is using a water pump as an impeller to grind 
air into water to create an air-water mixture that is sent to a porous hose to generate 
bubbles. The present invention is not limited to these two methods. It covers any 
methods of manipulating the local acoustical properties to construct a phase array. 
This would include, but is not limited to, air or gas, mixed or dissolved, in water or 
fluid, by mechanical, electrical, or chemical methods. A control (not shown) for 
activating the bubble diffuser 101 is attached to the source of compressed air. The 
control is used to activate the bubble diffuser 101 during the collection of seismic data 
using the sources 104 and receivers 105. 

As mentioned above, the bubble diffuser 101 is specially configured. 
Preferably, the bubble diffuser 101 comprises a plurality of diffuser elements 111. 
Embodiments of the diffuser elements 111 include, but are not limited to, ceramic 
stone air diffusers and porous hose air diffiisers. The diffuser elements 111 are 
preferably positioned in a vertically oriented phase grating. The diffuser elements 111 
are preferably positioned in the phase grating so that the emitted bubbles 108 
substantially suppress specular reflection of up-going seismic waves from the sources 
104. The phase grating is preferably focussed on a range of distances to the sources 
104 and receivers 105. 

The theoretical concept behind a first embodiment of the phase grating in the 
present invention is mathematically based on the concept of integer sequences from 
number theory. Hie mathematical basis for this is discussed in more detail below. 
These integer sequences can be used to construct the phase grating of diffuser 
elements 1 1 1 that will be used to interfere with multiples. Theoretically, the resulting 
phase grating diffuses the up-going seismic wave such that the multiples in the down- 
going specular direction are greatly reduced. For use in the present invention, the 
applicable types of integer sequences include, but are not limited to, quadratic residue, 
primitive polynomial, primitive root, complex Legendre, and Zech logarithm. In this 
first embodiment, the diffuser elements 111 of the phase grating are positioned 
according to a discrete number sequence. An example of this phase grating 
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configuration is shown in FIG. 2d. The preferred discrete number sequence is a 
primitive root sequence. The present invention also includes any type of phase 
gra ting other than those based on integer sequences or any other aspect of number 
theory. 

For instance, in another embodiment, the diffuser elements 111 of the phase 
grating are positioned in a vertically oriented, generally concave upward shape. An 
example of this configuration is shown in FIG. 2c. The preferred concave upward 
shape is a wedge shape, pointed upward. This wedge shape may be thought of as a 
continuous version of the discrete number sequence phase array described above. 

For the wedge-shaped embodiment of the bubble diffuser 101, the cross- 
section shape is maintained by the tension between the towing or auxiliary work 
vessels 106, 1 10 and buoys/stabilizeis 1 12. A stabilizer or clump weight (not shown) 
may be required at the base of the wedge to maintain the correct shape. For the 
primitive root embodiment of the bubble diffuser 101, a net, as shown in FIG. 2d, may 
be used as a support 107 to maintain the proper depths of the difiuser elements 111. 
The net is towed by towing lines 109 attached to towing or auxiliary work vessels 
106, 110 and buoys/stabilizers 112. Alternatively, for either embodiment of the 
bubble diffuser 101, the diffuser hose assembly may be attached directly to the towing 
or auxiliary work vessels 106, 1 10, rather than to the buoys/stabilizers 112. Since the 
bubble diffuser 101 is made of flexible components, such as cables, hoses, and net, the 
components can be stored on drums and easily deployed and recovered by the towing 
or auxiliary work vessels 106, 110. 

In order to achieve high acoustic reflection in the seismic frequency band with 
air, it is necessary to generate a layer of bubbles 108 (or air-water mixture) several 
meters thick A single line of diffuser elements 111 (ceramic stone, porous hose, or 
other) may not have a high enough airflow rate to achieve this. It is possible to 
increase the airflow rate by stacking multiple diffuser elements 111. FIG. 3a shows 
another embodiment of the wedge shaped bubble diffuser 101 with a plurality of 
bundles of porous hoses 301 stacked vertically. Although eight bundles of porous 
hoses 301 are shown in FIG. 3a, the present invention is not limited to that particular 
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number. To achieve a higher flow rate, a plurality of porous hoses 302 are wrapped 
around a supporting steel cable 303 in the circumferential direction as shown in FIG. 
3b. Although ten porous hoses 302 are shown in FIG. 3b, the present invention is not 
limited to that particular number. By a combination of concepts of the embodiments 
shown in FIGS. 3a and 3b, it is possible to increase the flow rate by fector of 8 x 10 = 
80. Another embodiment designed to increase the airflow rate is shown in FIG. 3c. 
In this embodiment, there are a plurality of trailing air lines 304 attached to a wedge 
shaped bubble diffuser 101. Although 300 trailing air lines 304 are proposed for the 
embodiment shown in FIG. 3c, the present invention is not limited to that particular 
number. Each trailing air line 304 emits bubbles 108 (as shown in FIG, 1). With 
horizontal extension in the inline dimension, the trailing air lines 304 can together 
feed more air than a single air line bubble diffuser 101. Since the trailing air lines 304 
are streamlined, this embodiment has less drag than the embodiment of FIG. 3a. The 
mflviTniim flow rate of this embodiment is easily adjusted by increasing the total 
length of the trailing air lines 304. As a variation of the trailing air lines, a continuous 
diffusing surface can be created by putting the trailing air lines adjacent to each other. 
Of course, all three of the concepts involved in FIGS. 3a - 3c can be combined to 
potentially increase the airflow rate by many orders of magnitude if necessary. 

FIG. 4 is a flow chart illustrating the steps of the method for one embodiment 
of the present invention for suppressing multiples in the collection of marine seismic 
data, using a specially configured bubble diffuser. The collection of marine seismic 
data takes place in a body of water with an air-water interface. First, in step 401, at 
least one seismic source is positioned in the body of water, as is typical in the 
industry. The seismic sources are preferably airguns or marine vibrators, positioned 
singly or in arrays, but may be any other appropriate marine seismic source known in 
the art. 

In step 402, at least one seismic receiver is positioned in the body of water, in 
the general vicinity of the seismic sources, as is typical in the industry. The seismic 
receivers are positioned below the air-water interface. The receivers are preferably 
pressure sensors such as hydrophones, but may be any other appropriate sensors 
adapted for the marine environment. 
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In a preferred embodiment, both the seismic sources and seismic receivers 
from steps 401 and 402, respectively, are towed through the body of water by seismic 
vessels equipped for this task, as is well known in the art The sources and receivers 
may be towed by the same or different vessels. The sources may be towed in single or 
multiple arrays and the receivers may be towed in single or multiple seismic 
streamers. In an alternative embodiment, the receivers may be positioned on the 
bottom of the body of water in an ocean bottom cable. Here the receivers are 
preferably pressure and particle velocity sensors such as hydrophones and geophones, 
respectively, but may be any other appropriate sensors adopted for the marine 
environment 

In step 403, a bubble diffuser is positioned in the body of water. The bubble 
diffuser is positioned so that the bubbles emitted from the bubble difiuser are 
positioned between the receivers and the air-water interface. The bubble difiuser is 
preferably towed through the body of water between the receivers and the air-water 
interface to properly position the emitted bubbles. The bubble diffuser may be towed 
by the same towing vessels towing the sources or receivers. Alternatively, the bubble 
diffuser may be towed by auxiliary work vessels, such as one or more tugs. 

In step 404, the bubble diffuser from step 403 is configured so that the emitted 
bubbles substantially suppress specular reflection of up-going seismic waves from the 
seismic sources. The details of this special configuration are described in conjunction 
with the discussion of FIGS. 2c and 2d above. In addition, the emitted bubbles are 
preferably of a sufficiently high volume concentration in the water to act as a high 
acoustic reflector of seismic waves. The emitted bubbles are also preferably of 
sufficiently small size that the bubbles rise slowly in the water, thus retaining their 
effective position for a longer time. 

In step 405, the bubble diffuser from step 403 is activated during the collection 
of marine seismic data by the sources and receivers. Thus seismic waves originating 
from the seismic source and traveling upward after possible reflection from reflectors 
below will be substantially deflected to the side rather than reflected down from the 
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air water interface toward the receivers. This will reduce the amount of noise caused 
by multiples as recorded by the receivers during the collection of marine seismic data. 

Mathematical Foundation 

The present invention depends upon the reflective properties of a phase array 
based upon an integer sequence from number theory. Among the possible phase 
gratings based on integer sequences, the primitive root sequence is known to 
minimize the main lobe and have evenly distributed side lobes. This is discussed in 
Schroeder, M. R., "Phase gratings with suppressed specular reflection", Acoustica, 
v.81, pp.364-369, 1995. Thus, the primitive root sequence is chosen to demonstrate 
how a phase array can reduce the multiples. 

In number theory, one can find a primitive root for any given prime number. 
The primitive root, g, for a prime number, p, is defined to be the smallest prime 
number g that generates a complete set of integers between 1 and p - 1 by remainders 
by the prime number p on the positive integer powers of g . Equation (1) shows the 
generation of a primitive root sequence. 

}nod( g J ,p),j- 1,2,... p Jnod( g 1 ,p),mod( g 2 lP )„..,mod( g p ~\p)}. (1) 

Here, mod(a t b) means the smallest positive remainder after integer multiples of b 
are subtracted from a. Thus, mod (6, 8) = 6-(0-8) = 6 and 
mod (1 1, 8) = 1 1-(1 • 8) = 3 . For example, prime number p = 1 has the primitive 
root g = 2 and together they generate the following sequence: 

{nod(2 1 1 11) l mod(2 2 > 11) mod(2 10 t 1l)}= {2,4,8,5,10,9,7,3,6,1}. (2) 

Using the primitive root sequence, a phase grating is constructed, with phases 
at different positions along the grating given by: 



♦,-^^,wh ere j=1A...,P-1. 



(3) 
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With the above phase grating, an expression of the linear array response A k in the far 
field, is obtained: 

A k =2exp(i(|) n+1 )exp(-27iinwsina k A) > (4) 

where i = V— 1 , X, is the acoustic wavelength in the medium, ct k is the incident polar 
angle (i.e. <x k = 0 corresponds to the normal or broadside direction), and w is the 
spacing between adjacent elements. By substituting the nondimensional spatial 

wavenumber k , defined as W ^ — sin a k , for the sine of the incident polar angle, 

Equation (4) becomes: 

A k =g exp(i^ 1 )exp[-2nink/(p-l)]. (5) 

The above expression is the discrete Fourier transform (DVT) of the phase grating 
contribution expOfy). 

The magnitude of ^ k is known when the nondimensional spatial wavenumber 
k is an integer: 



N- 



' 1 fork = 0 ^ 
-y/p for other integer values of k ' 



This expression simply states that the specular reflection (k = 0 ) is smaller than the 
side lobes (k * 0) by a factor of ^/p • Also, compared to the case of zero phase 
grating (i.e. perfectly reflecting surface with <|>j = 0 ), the specular reflection is 
reduced by a factor of p . Thus, when a large prime number is chosen for a primitive 
root phase array, one can greatly reduce the specular reflection and divert energy to 
the side directions. 
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Even though the primitive root sequence is shown as an example, other types 
of integer sequences or even real number sequences can be used in the present 
invention for multiple suppression. 

One way to create a phase grating for seismic (acoustic) waves is to vary the 
elevation of a reflecting surface. If the reference elevation of the reflecting surface 
(i.e. air-water interface) is set to be zero degree phase and the two-way path length 
between source and receiver is considered, the surface elevation Zj is related to the 
phase grating as: 



In Equation (7), X equals the wavelength for the lowest frequency in the signal band 
because the array cannot cover full angles below that frequency. 

FIG. 5 illustrates an example of a vertically oriented phase grating constructed 
in accordance with the present invention. Since the phase grating is discrete, the 
surface elevation is also discrete. Thus, the surface consists 'of flat reflecting surfaces 
501 with discrete elevation changes 502 and uniform spacing 503. The width 
(spacing) 503 of an individual surface can be found from the assumption behind 
Equation (4). In order for Equation (4) to be valid, the spacing w is: 



In Equation (8), X is the wavelength for the highest frequency in the signal band 
because it is required to use a more rigorous formulation than the array response, 
Equation (5), above that frequency. 

Equation (4) was derived assuming the distance from the array to a field point 
is much greater than the array dimensions. Thus, the phase grating is focused at 
infinite range and it may not perform well at finite ranges. In order to focus the phase 
grating at a finite range, the exact phase at a finite range should match the 




(7) 




\/2 for broadside incidence 
X/A for grazing incidence 



(8) 
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approximated phase at infinite range. For the case of normal incidence, the corrected 
surface elevation z' is obtained: 



(9) 



where r is horizontal coordinate at the array surface, R s is source range, and R r is 
receiver range. If z and r are small compared to the ranges, Equation (9) can be 
approximated as: 

2 '« 2+ I1Bb±Rl. (10) 
2 2R e R r V 

The above focusing technique is independent of phase grating type. 

In order to interfere with the multiple at the air-water interface, the size of a 
phase array must be large enough to interact with the Fresnel zone of the up-going 
seismic wave. The diameter of the Fresnel zone (FZ) is: 

FZ = ftXR B . (11) 

For a given sequence of phases for a phase grating, the number of periods of the 
sequence can be simply increased until the phase array covers the entire Fresnel zone. 

Examples 

Table 1 gives the input specification for an example of an embodiment of the 
present invention. This example is a design concept for a phase array multiple 
suppresser for a high-resolution seismic scale. The phase array is designed to 
suppress water-bottom generated multiples, including source and possibly receiver- 
side pegleg multiples. The receivers are configured as a streamer deployed below the 
multiple suppresser. 



Sound speed (c, mis) 



1500 
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Lower frequency ru-) 


200 


Upper frequency (f 2 , Hz) 


700 


Water bottom depth {d , m) 


100 


Source wavelet 


Ricker pulse with peak 
frequency at = 350 Hz 



TABLE 1 



A phase array is computed using two periods of a primitive root sequence for 
p = 31 . The profiles of the phase array are obtained and shown in FIGS. 6a and 6b. 
FIG. 6a shows the profile computed for this primitive root sequence according to 
Equation (7), while FIG. 6b shows the profile recomputed with the focusing technique 
of Equation (9). 

FIGS. 7, 8, and 9 show the theoretical performance of this phase array when it is 
subjected to a plane wave Ricker pulse. FIG. 7 shows the theoretical normalized 
amplitude 701 of the transfer function of specular scattering. In FIG. 7, the amplitude 
of the array transfer function becomes less than 0.2 above the lower design frequency 
(i.e. 200Hz) and decays as the frequency increases. It means that this phase array can 
effectively suppress the multiple by more than a factor of 5 above the lower design 
frequency. FIG. 8 shows the comparison between the theoretical spectra of specular 
reflection for the input pulse without the suppresser 801 and the spectra of the input 
pulse with the suppresser 802. The multiple suppresser clearly reduces the multiple 
energy throughout the design band. FIG. 9 shows the theoretical acoustic response in 
the time domain of the multiple without the suppressor 901 and the response with the 
multiple suppresser 902. The peak suppression ratio is about 24.3 (= 28 dB) and the 
root mean square suppression ratio is about 12.8 (= 22 dB). 

Although the example shown above is designed for high-resolution wavelengths, 
it carries over to standard seismic wavelengths (10 - 100 Hz frequency band), with a 
change in scale. Seismic-scale streamer lengths are 3 - 5 km. For a 100 m water 
bottom, and a 150 m wavelength, good multiple suppression is obtained for a diffuser 
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350 m wide around each streamer (or 175 m outboard of the furthest cable in a multi- 
streamer experiment). For a triple-streamer experiment with 100 m streamer 
separation in the cross-line direction, this implies a total diffuser width of 550 m. 
Diffuser maximum depth is half of the longest wavelength, which for 10 - 100 Hz 
seismic data would be 75 m. Most of the dements in the diffuser, however, would be 
shorter than this depth. For a 500 m water bottom, diffuser maximum depth does not 
change but its width grows to 700 m for the above triple-streamer embodiment 

Small-scale field experiments were performed to measure the actual performance 
of the above embodiment of the present invention. The experimental site is a 60 m 
square-shaped freshwater pond and the Tnavirmirn depth of the pond is 7 m. A 
watergun was used as a seismic source and 4 channels of hydrophone data were 
recorded. The air-generation mechanism used between 24 and 32 micro-bubble 
ceramic stones made by Point Four Systems, Inc., Canada (Model No. MBD-600). 
This particular ceramic stone can generate air bubbles as small as 100 microns in 
diameter. Two embodiments of the multiple suppresser of the present invention were 
deployed. The first configuration, shown in FIG. 10, is a simple wedge (V) shaped 
diffuser. The second configuration, shown in FIG. 11 is a primitive-root diffuser 
based on prime number p = 17 with finite depth focusing using the methods of 
Equation (9). The experimental setting miniaturized a common marine seismic data 
acquisition approach with a seismic vessel, towed streamers (linear hydrophone 
array), and an air-gun array (seismic source). The entire system (multiple suppresser, 
watergun, 4-channel streamer, and global positioning system) was assembled into one 
piece and towed by an electric winch. Each micro-bubble diffuser was fed with 
compressed air. 

Seismic data were collected and processed with this configuration. Seismic 
images were obtained both without and with the multiple suppresser operating, that is, 
generating the micro-bubbles. FIG. 12a is the baseline seismic image taken with the 
wedge shaped multiple suppresser of FIG. 10 and no multiple suppression. In the 
figure, the solid rectangular window 1201 indicates primary reflection of seismic 
waves from the water bottom and the dashed rectangular window 1202 corresponds to 
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the multiple energy generated by the primary reflection and the air-water interface. 
Those two seismic signals are highly correlated and thus the image is massively 
contaminated by multiple energy. FIG. 12b is another seismic image taken with the 
wedge shaped multiple suppresser and 40 psi compressed air supply to the micro- 
bubble generators. Here, there is multiple suppression by the method of the present 
invention. Compared to FIG. 12a, the primary seismic signal in the solid rectangular 
window 1203 remains unchanged, but the multiple in the dashed rectangular window 
1204 energy is visibly reduced. Even after the multiple time window, this image 
contains less reverberant energy, which is a high order multiple (multiple generated 
by a multiple). 

FIG. 13a shows the baseline seismic image taken with the primitive root multiple 
suppresser of FIG. 11 and no multiple suppression. As in FIG. 12a, the primary 
reflection exists in the solid rectangular window 1301 and the multiple energy exists 
in the dashed rectangular window 1302. FIG. 13b is another seismic image taken with 
the primitive roots multiple suppresser and 40 psi compressed air supply to the micro- 
bubble generators. Here, there is multiple suppression by the method of the present 
invention. The multiple energy in the dashed rectangular window 1304 is greatly 
suppressed compared to the baseline image in FIG. 13a. Even the primary reflection 
in the solid rectangular window 1303 is changed. While the baseline images in FIG. 
13a contains the receiver-side ghost 1305 (i.e. time-delay caused by the distance 
between the receiver and the air-water interface), the multiple suppresser interferes 
with the receiver-side ghost 1306 in FIG. 13b. It is additionally beneficial for 
interpretation of marine seismic data not to have the receiver-side ghost in the data. 

The small-scale field experiments verify that the method of the present invention 
is effective in removing multiple energy in marine seismic data acquisition. In order 
to implement the present invention in full-scale seismic acquisition, fast, easy, and 
reliable methods of deployment and recovery of the bubble diffuser are required and 
provided by the present inventioa Since the width of the multiple suppresser for full 
seismic scale may be on the order of a kilometer, flexible devices made of tension 
cable and air-diffusing hoses are preferred over solid/rigid devices. As discussed 
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above, FIGS. 2a through 2d show two example deployment concepts of the wedge- 
shaped and primitive root multiple suppressers for full-scale seismic acquisition. In 
addition to the main seismic vessel, two auxiliary work vessels assist to deploy, tow, 
and recover the multiple suppressers. Since the devices are made of flexible cables, 
hoses, and net, they can be stored on drums and easily deployed and recovered by the 
work vessels. Alternative configurations in addition to the embodiments discussed 
here are also possible. For example, the multiple suppresser may be towed in front of 
the seismic vessel, rather than between the seismic source and streamer. 

The present invention includes any other methods that deploy and recover the 
multiple suppressers, such as using neutrally buoyant and streamlined hydrofoils with 
encased ceramic air diffusers. The present invention also includes any other cable 
configurations that allow bubble-making devices to be attached to the cable and 
deployed. 

It should be understood that the invention is not to be unduly limited to the 
foregoing which has been set forth for illustrative purposes. Various modifications 
and alternatives will be apparent to those skilled in the art without departing from the 
true scope of the invention, as defined in the 
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CLAIMS 

What is claimed is: 

1 . An apparatus for suppressing multiples in the collection of seismic data in a 
body of water with an air-water interface, comprising: 

at least one source positioned in the body of water; 

at least one receiver positioned in the body of water below the air-water 
interface and near the sources; 

a bubble diffuser positioned in the body of water so that the bubbles emitted 
by the bubble diffuser are positioned between the receivers and the air- 
water interface, wherein the emitted bubbles provide high acoustic 
reflection and substantially suppress specular reflection of seismic 
waves; and 

a control for activating the bubble diffuser during the collection of seismic 
data by the sources and receivers. 

2. The apparatus of claim 1, further comprising: 

a support positioned in the body of water below the air-water interface and 
near the receivers and attached to the bubble diffuser; and 

a towing line attached to the support for the bubble diffuser for positioning the 
emitted bubbles between the receivers and the air-water interface. 

3. The apparatus of claim 1, further comprising: 

a compressed air hose for supplying compressed air to the bubble diffuser 
sufficient to emit bubbles providing high acoustic reflection. 

4. The apparatus of claim 1, wherein the bubble diffuser comprises a plurality of 
diffuser elements. 

5. The apparatus of claim 4, wherein the diffuser elements are ceramic stone air 
diffusers. 
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6. The apparatus of claim 4, wherein the diffuser elements are porous hose air 
diffusers. 

7. The apparatus of claim 4, wherein the diffuser elements comprise a plurality of 
porous hose air diffusers wrapped around the support 

8. The apparatus of claim 4, wherein the plurality of diffuser elements are 
positioned in a vertical orientation to increase bubble emission. 

9. The apparatus of claim 4, wherein the plurality of diffuser elements are 
positioned in a horizontal orientation to increase bubble emission. 

10. The apparatus of claim 9, wherein the plurality of diffuser elements comprise 
trailing porous air hoses. 

11. The apparatus of claim 4, wherein the bubble diffuser comprises diffuser 
elements positioned in a vertically oriented phase grating. 

12. The apparatus of claim 1 1, wherein the phase grating is focussed on a range of 
distances to the sources and receivers. 

13. The apparatus of claim 1 1, wherein the diffuser elements of the phase grating 
are positioned according to a discrete number sequence. 

14. The apparatus of claim 13 wherein the discrete number sequence is a primitive 
root sequence. 

15. The apparatus of claim 11, wherein the diffuser elements of the phase grating 
are positioned in a generally concave upward shape. 

16. The apparatus of claim 15, wherein the diffuser elements of the phase grating 
are positioned substantially in a wedge shape. 

17. A method for suppressing multiples in the collection of seismic data in a body 
of water with an air-water interface, comprising the steps of: 
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positioning at least one source and at least one receiver in the body of water, 
wherein the receivers are positioned below the air-water interface; 

positioning a bubble diffuser in the body of water so that the bubbles emitted 
by the bubble diffuser are positioned between the receivers and the air- 
water interface, wherein the emitted bubbles provide high acoustic 
reflection and substantially suppress specular reflection of seismic 
waves; and 

activating the bubble diffuser during the collection of seismic data by the 
sources and receivers. 

1 8. The method of claim 17, wherein the emitted bubbles are of sufficiently small 
size so that the bubbles rise slowly in the body of water. 

19. The apparatus of claim 17, wherein the diffuser elements comprise a plurality 
of porous hose air diffusers wrapped around the support. 

20. The apparatus of claim 17, wherein the plurality of diffuser elements are 
positioned in a vertical orientation to increase bubble emission. 

21. The apparatus of claim 17, wherein the plurality of diffuser elements are 
positioned in a horizontal orientation to increase bubble emission. 

22. The apparatus of claim 21, wherein the plurality of diffuser elements comprise 
trailing porous air hoses. 

23. The method of claim 17, wherein the bubble diffuser comprises diffuser 
elements positioned in a vertically oriented phase grating. 

24. The method of claim 23, wherein the phase grating is focussed on a range of 
distances to the sources and receivers. 

25. The method of claim 23, wherein the diffuser elements of the phase grating are 
positioned according to a discrete number sequence. 
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26. The method of claim 25 wherein the discrete number sequence is a primitive 
root sequence. 

27. The method of claim 23, wherein the difluser elements of the phase grating are 
positioned in a generally concave upward shape. 

28. The method of claim 27, wherein the diffuser elements of the phase grating are 
positioned substantially in a wedge shape. 
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Transfer Function of Specular Scattering : 1-D PR Grating 
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